Two impurities in a d-wave superconductor: local density of states 
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We study the problem of two local potential scatterers in a d-wave superconductor, and show 
how quasiparticle bound state wave functions interfere. Each single-impurity electron and hole 
resonance energy is in general split in the presence of a second impurity into two, corresponding 
to one even parity and one odd parity state. We calculate the local density of states (LDOS), and 
argue that scanning tunneling microscopy (STM) measurements of 2-impurity configurations should 
provide more robust information about the superconducting state than 1-impurity LDOS patterns, 
and question whether truly isolated impurities can ever be observed. In some configurations highly 
localized, long-lived states are predicted. We discuss the effects of realistic band structures, and 
how 2-impurity STM measurements could help distinguish between current explanations of LDOS 
impurity spectra in the BSCCO-2212 system. 



I. INTRODUCTION 

The study of isolated point-like impurities in li-wave 
superconductors began with the observation that the 
scattered wavepattern probed in tunneling experiments 
should inherit the fourfold symmetry of the d-wavc 
stateQ. It was pointed out shortly afterwards by 
Salkola et aZ.0, following earlier work on the analo- 
gous p-wave problem by StampQ, that quasi-bound res- 
onances should occur near strongly scattering impuri- 
ties. Both the resonance energy and the details of the 
spatial structure of the resonant wave-pattern depend 
sensitively on electron correlations and a large body 
of theoretical work exists exploring this relationship jj]. 
Assuming that the important details of the scattering 
potential are understood, careful studies of the local 
density of states (LDOS) near isolated impurites pro- 
vide a uniquely powerful probe of the superconducting 
state. Experimentally, the high temperature supercon- 
ductors (HTSC), notably Bi 2 Sr 2 CaCu 2 8+I 5 (BSCCO), 
have been examined with sensitive scanning tunneling 
microscopy (STM) techniques. Images of the local envi- 
ronment of impurities [5|, |fj have confirmed the existence 
of quasi-bound states near strongly scattering impurities 
like Zn, but have led to new questions regarding the mi- 
croscopic model for impurities and the superconducting 
state itself. 

Other types of inhomogeneities not directly correlated 
with impurity resonances have been discovered by STM 
measurements on the cuprates. These include the ob- 
servation of large, quasi-bimodal nanoscalc fluctuations 
of the superconducting order parameter, together with 
spatially correlated variations of the electronic struc- 
ture and lifetime. 0, @ More recently, the observation 
of checkerboard patterns in the LDOS of vortex cores[^ 
and in inhomogeneous samples in zero fieldB3 in BSCCO 
has led to speculation that antiferromagnetic phases 
can be formed in regions where superconductivity is 
supressed. ^lj A good deal of subsequent theoretical work 



on the competition between ci-wave superconductivity 
and various exotic order parameters nas excited the 
high-T c community with the suggestion that STM mea- 
surements could be revealing the presence of competing 
magnetic or other exotic subdominant order, shedding 
light on the origins of superconductivity itself. 

On the other hand, a more conventional but still fasci- 
natin g ex planation has been put forward by by Hoffman 
et al.[13, who argue that Friedel oscillations related to 
quasiparticle scattering between nearly parallel sections 
of the Fermi surface can lead to the observed checker- 
board patterns, which correspond to well-defined dynam- 
ical peaks in momentum space. These general remarks 
were followed very recently by a calculation by Wang and 
Lee0 studying the Friedel oscillations of a single im- 
purity in a d wave superconductor in momentum space. 
These calculations support the idea that what STM is 
seeing is primarily the effect of quasiparticle wavefunc- 
tions interfering with one another in a fluctuating disor- 
der field on a <i-wave superconducting background. 

Studies of true interference of quasiparticle wave func- 
tions in the presence of more than a single impurity are 
rare, however. Numerical solutions of the many-impurity 
problem yield little insight into the mechanism of inter- 
ference itself. In principle, the problem of two impurities 
is the simplest one which displays the interference effects 
of interest for the present problem; we study it here in 
order to understand whether the STM analysis of iso- 
lated impurity resonances is indeed justified, and to test 
if the scenario proposed by Hoffman et al.[L| is tenable 
in the presence of interfering Friedel oscillations. There 
is a small amount of earlier work on two impurities in a d- 
wave superconductor, most of it also numerical. Onishi et 
al.[l5j solved the Bogoliubov-de Gennes (BdG) equations 
and presented a few local density of states profiles for two 
impurities. In a work philosophically related to ours, the 
interference of bound state wave functions in the fully dis- 
ordered system was discussed by Balatsky in the context 
of supression of localization effects due to impurity band 
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formation[16j. Micheluchi and Kampf have recently ex- 
hibited numerically how impurity induced bound states 
accumulate at low energies, and argued that the im- 
purity band at low energies could be studied from this 
perspective. Finally, during the preparation of this 
manuscript, a paper by Morr and Stavropoulos exam- 
ined the two-impurity problem with particular emphasis 
on predictions for cuprate STM studies. [3 We comment 
throughout the text on comparisons with these previous 
works. 

We begin in Section ^ by reviewing the solution to 
the single impurity problem which has been studied by 
several authors. In Section IlIII we set up the formalism 
for the two-impurity problem and give the exact solution 
for the T-matrix, as well as a simplified form for some 
special cases. In Sec. IIVI we discuss the dependence 
of the resonance splittings on the orientation and magni- 
tude of the intcrimpurity separation R. The resonant en- 
ergy splittings are argued to give important information 
which is qualitatively different from that obtainable from 
1-impurity resonance energies; in particular, they allow 
one in principle to map out, by repeated measurements 
of energies for impurity configurations at separation R, 
the spatial structure of the homogeneous superconduct- 
ing Green's function. Starting from two widely separated 
impurities with R ^> £o, where £o is the coherence length, 
we then give analytical solutions for the splittings which 
exhibit explicitly the strong dependence on the direction 
of R due to the d-wave nodes. Impurities located along 
a 100 axis interact very weakly for R > t;o, whereas con- 
figurations near 110 relative orientation lead to strong 
hybridization. 

In Section we begin by discussing the simple quan- 
tum mechanics problem of how bound eigenstates with 
two scattering centers are constructed from the eigen- 
states of one, and show how these states may be classi- 
fied. Of four states which arise from the single impurity 
particle and hole resonances, two are spatially symmetric 
(s), and two antisymmetric (p), one on each side of the 
Fermi level. The energy ordering of these states, as well 
as the crude qualitative interference pattern (construc- 
tive or destructive) depends on the configuration R in a 
predictable way. The local density of states (LDOS) is 
then calculated and plotted for several cases to illustrate 
the spatial dependence of the resonant state wave func- 
tions. Spectra on individual sites reveal unexpected phe- 
nomena: in certain circumstances the 1-impurity states, 
which become sharper as they approach the Fermi level 
due to the coupling to the linear e?-wave continuum, inter- 
fere to create localized, extremely sharp states located at 
energies quite far from the Fermi level. In Section IVll we 
discuss the situation for a more realistic band character- 
istic of the BSCCO-2212 system on which impurity STM 
studies have been performed. The trapped quasiparticlc 
states are still found, and a dynamical rcsonanance crite- 
rion depending on the exact Fermi surface and impurity 
orientation R is identified. 

Finally, in section VII we present our conclusions and 



discuss the implications for STM experiments and other 
aspects of the disordered quasiparticlc problem. 



II. SINGLE IMPURITY PROBLEM 

In broad terms, there are two distinct points of view 
which have evolved regarding the localized scattering res- 
onances observed by STM on the surface of BSCCO-2212. 
In the traditional quasiparticlc picture, one simply cal- 
culates the scattering T-matrix for non-interacting BCS 
quasiparticles in a d-wave superconductor and finds a 
pair of resonances which are approximately symmetric 
in energy about the Fermi level, and which have distinc- 
tive spatial patterns^]. In the second point of view, ad- 
ditional physics arises from the local disruption of the 
strongly-correlated ground state by the impurity which 
is predicted to break singlet correlatio ns [Hi l20| . nucleate 
short-ranged antiferromagnetic order|2ll l22l l23l 123. 12M 
l2t| . or spin fluctuations associated with a nearby phase 
transition [Til |27| . Perhaps the strongest motivation for 
this point of view is the observation of local moments 
near Li and Zn impurities in the superconducting state 
of underdoped YBCO, as observed in NMR,[2|, ^ ex- 
periments. For the STM experiments, the basic ques- 
tion is whether the spatial structure which is observed 
at low energies is effectively the result of BCS quasipar- 
ticles scattering from a short-range potential (ie. Friedel 
oscillations), or requires consideration of local correlation 
effects (eg. spin density wave formation|27j) or dynamical 
effects (Kondo physics |l9ll30j). 

One notable feature of the STM experiments on Zn- 
doped samples of BSCCO is that only a single, nega- 
tive energy resonance with a resonance energy of f» —1.5 
mcV is seen;0 the predicted impurity-induced reso- 
nances come in mirror pairs as a consequence of the 
particle-hole symmetry of the superconducting state. A 
further inconsistency is the fact that a large LDOS is ob- 
served on the Zn impurity site; the strong repulsive po- 
tential which Zn is believed to possess must necessarily 
allow little or no electron spectral weight at the impurity 
site. One appealing explanation |3 ll l32l | is that the mea- 
sured LDOS of the Cu02 planes is in fact filteredby an in- 
ert surface layer, leading to an apparent redistribution of 
spectral weight in the tunneling LDOS. With this mech- 
anism, one can simply understand the LDOS without in- 
troducing strong correlation physics. We note, however, 
that while this mechanism explains the observed LDOS 
for Zn impurities, it is problematic for both Ni impurities 
and Cu vacancies, which are consistent with the quasi- 
particlc picture without invoking a filtering mechanism. 
(To date, there is no convincing model which has ex- 
plained the spatial distribution of the LDOS in all three 
cases.) One of the goals of this work is to study the effect 
of the filtering mechanism on the resonant structure of 
two closely-spaced strongly-scattering impurities within 
the quasiparticlc point of view, providing a more rigor- 
ous test of the quasiparticlc-plus-filtcr mechanism for Zn 
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impurities. 

The BCS Hamiltonian for a pure singlet superconduc- 
tor can be written as: 



(1) 



where $k = (ckj,c^ k ^), is a Nambu spinor. In this 
work we will consider several forms for the dispersion 
ek- Analytic results are presented for a parabolic band 
6k = k 2 /2m, with corresponding <i-wave order param- 
eter Ak = A maa: cos 2(f> ((f> is the angle in momentum 
space which k makes with the 100 axis). Numerical 
results are presented for a simple tight binding model 
£k = — 2t(cosk x + cosky) — fi and for a realistic 6- 
parameter tight-binding model proposed by Norman et 
al. |33j| . both having the corresponding c?-wavc order pa- 
rameter Ak = Ao(cosfc x — cosky). Note the maximum 
value of the order parameter in the lattice system with 
the current convention is 2Ao. The matrices Tj are the 
Pauli matrices. 

The Hamiltonian of a single on-site impurity at r = 

may be written as Hi mp = Ek k' Mj^k^^k' i w here 
Vo is the strength of the impurity potential. The 
Green's function G kk ' ( w ) m the presence of the impu- 
rity is expressed in terms of the Green's function G k (aj) 
for the pure system as G(k, k ,w) = G°(k, w)<S kk ' + 
G°(k, w)T(w)G°(k ,oj), where the " symbol indicates a 
matrix in Nambu space. The solution is 



T = TqTq + T 3 T 3 

To = V 2 G /(S+S-) 

T 3 = V 2 (c-G 3 )/(S + S-), 



(2) 



where Go and G 3 are the To and t 3 Nambu components of 
the integrated bare Green's function G°(k, u>). This 
expression has resonances when 



S± = 1-V (G 3 tG ) = 0. 



(3) 



Note that in a simple band l/(7rAr Vo) = c/(irN ) is the 
cotangent of the s-wave scattering phase shift r/o, where 
iVo is the density of states at the Fermi level. In the spe- 
cial case of a particle- hole symmetric system, G3=0 and 
the resonance energy is determined entirely by Go, which 
is given in the case of a circular Fermi surface by Gq(u>) = 



2jt 



U)\U) 



A k ] 1 I 2 which for low energies u> <C A, 



takes the form G (uj) ~ -(■Ku/A max )(log4A max /u + i). 
One may then solve Re S±(w + i0 + ) = and estimate 
the resonance width T on the real axis. In the case of 
strong scattering c <C No, the resonance energy and 
scattering rate T are 



r = 



±7rcAo 
2 log(8/7rc) 

7T 2 CA 

41og 2 (8/7rc) 



(4a) 
(4b) 



This result was first obtained by Balatsky et al.Q, fol- 
lowing earlier work on the p-wavc analog problem by 
Stamp j3] . Note that the resonance becomes a true bound 
state only exactly at the Fermi level £1 = 0, when c = 0; 
for finite c there are two resonances whose energies are 
symmetric, fig" = —^Iq m this approximation. As seen 
from ©, in particle-hole asymmetric systems (G 3 7^ 0), 
the resonance is tuned to sit at the Fermi level for some 
value of the impurity potential Vo which is not infinite, so 
the term "unitarity" (as used in this work, f2 = 0) and 
"strong potential" (Vo — > 00) are not synonymous. |34l3a | 

The ambiguity in defining the resonance energy pre- 
cisely arises already at the level of the 1-impurity prob- 
lem. Eq. @ is in fact an equation for a complex fre- 
quency u>, which may be shown to have no solution in 
the upper half plane. Thus the T-matrix has no true 
pole at any u) = O' + ifl" in the complex plane, but only 
a maximum which lies along the real axis. Only when the 
real part f2' approaches the Fermi level does the damp- 
ing become sufficiently small to allow one to speak of a 
well-defined resonance, in this case the real part f2' ap- 
proaches the solution Oq of Re[S±] = for uj on the 
real axis. For this reason the "resonance energy" is usu- 
ally taken to be fig , as given, e.g. in Eq. I|4a|) for the 
particle-hole symmetric case. It is important to keep in 
mind, however, that the definition becomes meaningless 
as the bound state energy moves far from the Fermi level; 
for example, the apparent divergence of l|4a(l as c — > 8/ir 
is artificial, since no well-defined resonant state exists by 
the time f2o is a significant fraction of the gap Ao. It 
should also be noted that the generalization of this res- 
onance criterion to more complicated situations, where 
the denominator does not factor, is not straightforward. 

Even if the energies of the particle and hole resonant 
states are symmetric, their spectral weights on a given site 
may be quite different. 0jQ The finite impurity potential 
acts as a local breaker of particle-hole symmetry, leading 
in the case of repulsive potential to a large peak in p(r, u>) 
at the impurity site r = at negative energy (holclikc 
states) and a small feature at positive energy (electron- 
like), as seen in Fig. ^ The impurity-induced LDOS 
decays as r~ 2 along the nodal directions (for the particle- 
hole symmetric system), and exponentially along the 
antinodes. The LDOS in the near field is more compli- 
cated, however: the nearest neighbour sites have peaks 
at ±f2o, with the larger spectral weight at +^0- In the 
crossover regime r ~ £o> the LDOS is enhanced along 
the node direction for holelike states, but is spread per- 
pendicular to the node direction for electronlike states. 
These spatially extended LDOS patterns are the finger- 
print of the impurity-induced virtual bound states. In 
Figure ^ we illustrate the LDOS pattern expected for 
both particles and holes for a resonant state close to the 
Fermi level. Results are obtained with a simple half-filled 
tight-binding band, £k = — 2t(cos k x +cos k y ), and unless 
otherwise specified all energies are given in units of the 
hopping t. 

At the present writing, the LDOS pattern produced 
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FIG. 1: Summary of LDOS results for 1-impurity problem 
on a tight-binding lattice, Ao = 0.1, /x = 0, Vo = 10, |Qq | ~ 
0.013: a) LDOS vs. uj on the impurity site; b) LDOS vs. uj on 
nearest neighbor site; c),d) LDOS map at resonance uj = Qq 
and uj = Color scales in c) and d) are relative to the 

nearest neighbor peak heights shown in b). 



within the simple T-matrix theory for a single strong im- 
purity, while 4-fold in symmetry, does not agree in detail 
with STM experiments on Zn impurities and native pla- 
nar defects. [3 It is not currently clear whether this is due 
to a failing of the microscopic impurity model, e.g. fail- 
ure to include strong correlations in the host or magnetic 
degrees of freedom or whether a relatively trivial tunnel- 
ing matrix element effect prevents direct observation of 
the simple pattern by STM [3iTI3^ . For the moment then 
we consider only the simplest 2-impurity model possible, 
recognizing that direct application to experiments awaits 
a resolution of the discrepancy at the 1-impurity level. 



III. T-MATRIX FOR 2 IMPURITIES 

We now introduce the formalism necessary to study 
the interference between two resonances of the type dis- 
cussed above when two identical impurities are brought 
close to one another. The perturbation due to 2 identical 
impurities located at positions R^ and R ?n , is given by 



Hi, 



To 



E 



(5) 



where 4?^ = [c^qJ. By iterating the procedure for the 
single impurity T-matrix with two impurities present we 
find, in a 4 x 4 basis of spin and impurity site labels: [3l| 



ff m G°(-R,)f e 



ff e G°(R)f n 



(6) 



where R = R^ — R m and where Tg , T m are the single 
impurity T-matrices associated with the two impurities. 
For identical impurities, Tg = T m = T(ui), the single 
impurity T-matrix defined previously. The quantity / is 
defined as: 



f(oj) = [1 - G (-R, W )f,MG°(R, w )T m M]- 



(7) 



where G°(R, uj) = ^ k exp[ik-R]G^(o;) is just the Fourier 
transformation of G k (tj), the unperturbed Nambu 
Green's function. For systems with inversion symmetry 
G°(R,w) = G°(-R,w). Note that in Eq. ©, the phys- 
ical processes are clearly identifiable as multiple scatter- 
ings from each impurity £ and m individually, plus in- 
terference terms where electrons scatter many times be- 
tween I and m. In k-space, we can write the T-matrix in 
the more usual 2x2 notation as 



7WH = [e ik -^r e*'*™T ]T t 



e - 4 k'.R iTo 
e - 4 k'-R mTo 



(8) 



where To is the Pauli matrix. In Sec.^we showed that, 
provided the resonance energies are distinct, peaks in the 
total density of states correspond to minima of the T- 
matrix denominator: 

V = det[l - G°(-R, w)f (w)G°(R, w)f(w)]. (9) 

Explicitly, V = V{D 2 I{S\S 2 _) with 



V 2 = V+Vz 



V 2 GI(R,uj) 
V 2 G 2 (R,uj) 



(10) 



where 



V, 



[l-V G 3 (0,w)±VoGo(0,w)] 
+(-l) a K)[TG (R, W ) + G 3 (R,^)] 



(11) 



The factors T>\ , T> 2 determine the four 2-impurity res- 
onant energies. Here G Q (R, uj) is the r Q component of 
the integrated bare Green's function 



G Q (R, uj) = ^Tr{T a G°(R,uj)) 



(12) 



For completeness, we also give the explicit form of the 
T-matrix itself: 



2"k,k' (w) 

where 

M P -- 
M, -- 



2V 



cos(k • R/2) cos(k' • R/2) M s (uj) 

(13) 



+ sin(k • R/2) sin(k' • R/2) M p (oj 



V+V x G 1 (R,oj)V V 1 
G^^VoVt -V 2 ~V 1 

V+V 2 -G 1 (R,uj)V V 2 
-G x {R,uj)V q V 2 -V^V 2 



(14) 
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In certain special configurations, e.g. if the two impuri- 
ties are located at 45° with respect to one other, it is easy 
to check that Gi (R, ui) — V R. In this case the entire 
resonant denominator factorizes T> = T>i + T>i-T>2+'D2- ■ 
The T-matrix then takes the simple diagonal form 



Tk, k >H = 2y cos(k-^)cos(k'-5:) 



+ 2V sin(k • 5l) sin(k' .5) 



T+ 




Vi- 








V 2 - 





(15) 



where t± = (T3 ± tq)/2. 



IV. BOUND STATE ENERGIES FOR TWO 
IMPURITIES 

Measuring bound state energies of impurity resonances 
in STM experiments allows one to obtain information on 
impurity potentials, and has the virtue of being indepen- 
dent of the STM tunnelling matrix elements. On the 
other hand, resonance energies of isolated single impu- 
rities provide no information on the spatial structure of 
resonant or extended state electronic wavefunctions. In 
principle, measurement of only the resonance energies of 
isolated pairs of impurities with different separations R 
is the simplest method of getting spatially resolved infor- 
mation on electronic wave functions independent of the 
exact tunnelling mechanism. 

When two identical impurities with resonance energies 
£Iq , £Iq are brought together, the bound state wavefunc- 
tions interfere with one another, in general splitting and 
shifting each resonance, leading to four resonant frequen- 
cies fli , £lf , fl 2 and fij, where the subscript indicates 
which factor in Eq. Ijllfl is resonant. If splittings are not 
too large, the electron and hole resonances are related in 
a similar way as in the 1-impurity problem, fi^ ~ — 
and Vt 2 ~ — ^2 ■ Again the weight of each resonance 
may be quite different or even zero on any given site. A 
large splitting may be taken as evidence for strong hy- 
bridization of quasiparticle wavefunctions. If we take the 
interimpurity distance R as a parameter and keep impu- 
rity potentials and other parameters fixed, there are two 
obvious limits where this splitting vanishes. In the case 
of separation R = 0, the two impurities combine (math- 
ematically) to create a single impurity of strength 2Vb, 
so both flf 2 approach the f^(2Vb) appropriate for the 
double strength potential. In the case of infinite separa- 
tion R — ► oo, we must find f2^" 2 approaching the Oq(Vo) 
appropriate for isolated single impurities. 



A. Gas model 

Eq. © is a general result for two (5-function potentials 
embedded in a host described by an arbitrary Gq. We 



would like to derive analytical results for the resonance 
energies obtained therefrom to get some sense of the ap- 
propriate length scales and symmetries in the problem. 
At large distances, the resonance energies must approach 
the single impurity values, so the splittings can be cal- 
culated perturbatively. To do so one must first obtain 
analytical expressions for the large-distance behavior of 
the unperturbed Green's functions. This is difficult for 
the superconducting lattice tight-binding model on which 
most of this work is based, but much insight can be 
gained by studying the equivalent gas model, with spec- 
trum ek = k 2 /2m. In this case expressions have been ob- 
tained by Joynt|34| and Balatsky et al.^(| for the d-wave 
integrated Green's functions G Q (R, oj) at large distances, 
both for R making an angle 45° or 0° with the x axis. 
For frequencies w/A <C 1/kpr <C l/kp£o, these reduce 
to 



G°(R,cj) 



(16) 




N e 

JV ° k F R 4+7r 2 £gfc 



T"0 + Ti + T 3 ) COS UfR- 

ti — T3) sin kpR 



R 
R 



(110) 
(100) 



where £0 — vp /nA max is the coherence length. 

The resonance energies may now be found by inserting 
these expressions for frequencies to = f2^ + 5 into © and 
solving for the shifts S. We find £lf 2 — ^0 ^ ^ with 



(17) 



log(O /A m 



sin UfR fcf^o 



/feFRlog(f2o/A mM ) 



kF R A+^r R II (HO) 
cos(/c F i? + 7r/4) R II (100) 



These expressions are valid for 8/Q.q <C 1. 

Clearly the decay of the splitting ~ cxp —r/£o/y/kpr 
is much more rapid for distances larger than the coher- 
ence length along the antinode (100) than for along the 
nodes, where it falls as ~ 1/r. The lack of a scale in 
the long-distance interference of quasiparticle wavefunc- 
tions oriented along (110), where they strongly overlap, 
is of potentially crucial importance in the STM analysis 
of "isolated" impurities, and we will bear this question 
in mind in what follows. 



B. Lattice model 

Here we consider a tight-binding model for ease of nu- 
merical evaluation. The definition of the resonant ener- 
gies can be obtained either by finding the minimum of 
© or from an analysis of phase shifts [35j|. The solutions 
corresponding to each factor in Eq. (|10fl can then in gen- 
eral be tracked as a function of separation R by minimiz- 
ing 2?i t 2 separately. In practice, this works well except in 
some special cases where the minima are very shallow. In 
Figurc[5|wc show the result for a particle-hole symmetric 
system. It is seen that each factor T> a corresponds to 
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FIG. 2: 2-impurity resonance energies f2^ 2 > vs. impu- 
rity separation R/(y/2a) for R || (110), /i=0, A = 0.1, 
Vo = 10. Solid line: upper 2-impurity positive resonance 
energy; dashed line: lower 2-impurity resonance energy; up- 
per dashed-dotted line: reference 1-impurity resonance energy 
Qq (Vo)- Lower dashed-dotted line: 1-impurity resonance en- 
ergy for double impurity strength S7^~(2Vb). Symbol indicates 
resonant channel: o—Q,^ > 0; *=f2j > 0. 
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FIG. 3: 2-impurity resonance energies Qf 2 /t > vs. im- 
purity separation R/(a) for R || (100), /i=0, Ao = 0.1, 
Vo — 10. Solid line: upper 2-impurity positive resonance 
energy; dashed line: lower 2-impurity resonance energy; up- 
per dashed-dotted line: reference 1-impurity resonance energy 
f2g"(Vo). Lower dashed-dotted line: 1-impurity resonance en- 
ergy for double impurity strength S7j(2Vo). Symbol indicates 
resonant channel: o=f2+ > 0; -k—fl^ > 0. 



an oscillating function of R, with the factor determining, 
e.g. fij, changing from site to site according to whether 
the site is even or odd. This is due to the strong R de- 
pendence of the components G a ; in the simplest case, 
R || (110) and n = 0, G 3 (R,u) = Gi(R, J) = but 
G (R,w) = J2k cos(k x R/V2)cos{k y R/V2)Gl(Lu) oscil- 
lates rapidly. At R = 0, the problem reduces to the 
double-strength single impurity case; the factor 2?i gives 
the resonant frequency f2^(2Vb) and the factor T>2 is 
1. At large separation the tt^ and tt^ "envelopes" are 
seen to converge to fi^(Vb) with a length scale of a few 
£o — 10a for the parameters chosen. 

In the R || (100) case, the oscillations of the bound 
state energies with increasing R are not so simple, as 
seen in Fig. [3J The one obvious simple difference from 
the (110) case is that the energy splittings vanish much 
faster with distance, as expected from the discussion in 
Section TlV Al Otherwise the short distance behavior of 
the bound state energies is complicated. One can check 
that the energy closest to the Fermi level is fi^ when 
R = 2 + 4n, n integer, and otherwise. 

In general, the short distance behavior is difficult to 
analyze analytically and we note that in neither the 
(110) or (100) direction do resonances appear at all for 
R = 1. Clearly the hybridization is so strong in these 
cases that the picture of perturbatively split 1-impurity 
states breaks down. More importantly, the splittings are 
significant out to quite large distances. Parameters in 
Figures [5] and [3] are chosen such that £o f=a 10a, as seen 
from Figure 13 where we indeed expect a e~ R ^° falloff 
according to the previous section. On the other hand, 
Figure [2] indicates strong interference out to distances of 



30a or more! 



V. LOCAL DENSITY OF STATES 
A. General 

The spatial distribution of two-impurity bound states 
gives much more detailed information about the nature 
of the quantum interference processes between impuri- 
ties than the bound state energies by themselves. Here 
we ask whether one can simply express the bound state 
wave functions of the 2-impurity system in terms of the 
1-impurity bound states, and how they can be classi- 
fied by symmetry. We would like to make predictions for 
STM experiments, including which qualitative features of 
the spectra reflect the quantum numbers of these states 
directly, and how these features depend on impurity po- 
tential and configuration. As indicated in the previous 
section, it is important to determine how far apart two 
impurities need to be to be considered "isolated" . Fi- 
nally, we would like to understand how robust these pre- 
dictions are with respect to changes in band structure, 
scattering potential, etc. 

Throughout this work, the LDOS refers to the tunnel- 
ing density of states, 

^( r ' w ) = \ < YjPa{ T > w ) ( 18 ) 
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with the spin-resolved LDOS, 

Pt(r,w) = -7r -1 Im G u (r,r, w + i0 + ) (19a) 
^(r,^) = +7r" 1 Im G 22 (r,r,-w-i0 + ) (19b) 

where the subscripts 11 and 22 refer to the electron and 
hole parts of the diagonal Nambu Green's function. 

The factor of 1/2 in Eq. |JTHJ) ensures that the LDOS 
normalization is 



du>p(r, u>) = 1. 



Provided the peaks in the LDOS are well defined, the 
peak energies agree closely with the resonance energies 
defined by Eq. I|4a|) . Some care is required, however, 
because peaks in the LDOS may not appear on all sites, 
and can be difficult to resolve. It is also worth noting that 
the peaks in the LDOS are not symmetric with respect 
to the Fermi energy, though in practice the degree of 
asymmetry is very small. 

The local electron density of states measured by STM 
is given by Eq. (|18fl . with 

G(r,r,w) = ^G°(k,w) (20) 

k 

+ e- i;(k - k,) - r G°(k, W )2W MG°(k', «,), 

k,k' 

and Tkk' is the T-matrix for any number of impurities. 

B. Interference of 1-impurity wavefunctions. 

In the 1-impurity case, the T-matrix is given by Eq. 
(0) and it is easy to see that 



<JG&(r,r,w) = -V lm 



(G^r)) 2 , (G° 21 (r)f 



s + 



(21) 



Quite generally one can express the Green's function in 
terms of the exact eigenstates V'ri(r) of the system in the 
presence of the imnuritv|37| 



<5Gn(r,r, u) 



C(r)^»(r) 

n 

U! — fir. 



(22) 



where the final approximation is valid for a true bound 
state with u> very close to a particular bound state energy 
f2„, and will be a good approximation in the present 
case to the extent the resonances are well defined, in the 
sense discussed above. Comparing with the form 121|) 
thus allows us to identify the positive and negative energy 
wavefunctions of the single-impurity resonances (Vo > 
assumed): 



ip±(r) 



G§i(r,w) u>- 



G?i(r,w) w 



W 



(23) 



where Z± are non-resonant wave function normaliza- 
tion factors. Note that the electron-like bound state 
eigenfunction is directly related to the off-diagonal bare 
Green's function, while the hole-like wave function is pro- 
portional to the diagonal bare Green's function. 

We can follow the same procedure for the 2 impu- 
rity Green's function, and ask how the eigenfunctions 
at a particular resonant energy are related to the sin- 
gle impurity wave functions we have just found. Since 
the single-impurity resonant energies arc different from 
the 2-impurity energies, this analysis will be valid to the 
extent the splittings are small compared to f2o±- The 
Green's function <5G(r, r) can now be constructed from 
Eq. (|13|> an d the wave functions read off by comparing 
with the spectral representation in the same way as in 
the 1-impurity case. By examining (|13|) it may be shown 
that, depending on whether T>\ or T>2 is resonant, the 
wave functions thus extracted will be of definite spatial 
parity, ip n (r) = ±tp n (-r). We find 



z p _ (g\ 



r_ = z s _ ( g° u 



Gi(R)Vb ^ 
V+ 12p 



lp Gi(R)K, 2p 

ys f flO g l( R )K) ^,0 



2?r 



Gi(R)Vc 



LU = 0,2- 

u> = f2i+ 

ui = f2i_ (24) 



where G? > 

7S,P 



G°(r - R/2) ± G°(r + R/2), and the 
Z± p are normalization coefficients. These are the two- 
impurity odd (p) and even-parity (s) resonant state 
eigenfunctions expressed directly as linear combinations 
of the corresponding one-impurity eigenfunctions tJj± 
given in l|23|l . 



1. R || (110) 

We note now that in general particle and hole-like 
1-impurity eigenfunctions are mixed in each 2-impurity 
state (|24|) : this is possible because anomalous scattering 
processes with amplitude Gi(R) can take place. There 
are special situations, including all configurations with 
R || (110), where Gi(R) = and the eigenfunctions be- 
come much simpler and do not mix particle and hole 
degrees of freedom, 



V4M = Z% 



G? ls (r,u/) u = Sl 1 . 



(25) 



Re if Im¥ l^l 2 LDOS 




FIG. 4: Comparison of symmetric (s) and antisymmetric (p) 
combinations of 1-particle wave functions with exact LDOS 
for Ao = 0.1, Vo = 10, fx = 0, R = (6,6). 2-impurity reso- 
nance energies are fl2±/t = ±0.0195, Q.\±/t — ±0.0075. ener- 
gies are always ordered from highest (top) to lowest (bottom) . 




FIG. 5: Comparison of LDOS spectra p(r,u) for various r 
with two impurities with Ao = 0.1, Vb = 10, \i = at positions 
(—3, —3) and (3, 3) (R = (6, 6)). Vertical lines correspond to 
1-impurity resonances at fi^ = ±0.013. Note difference in 
vertical scales between upper and lower panels. 



and 

We study this case now both to get a simple idea of how 
the quantum interference between 2 impurities works in 
practice, and because the (110) states are claimed to be 
of particular importance for the formation of the impu- 
rity band due to the strong interference of long range 
tails in the quasiparticlc wavefunctions. In Figure 01 
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FIG. 6: LDOS maps at resonant energies for R || (110). 
Tight binding band, A = 0.1, V = 10, (J, = 



we have exhibited the resonant state wavefunctions Tp± p 
for the (110) case. We exhibit both the wave functions 
themselves, so that the reader may verify the parity of 
these states explicitly, and their absolute square. There 
is good agreement between the spatial pattern of the \ip\ 2 
as defined and the exact LDOS calculated from i|21l ) at 
each resonant energy, implying that near each resonant 
energy the nonresonant contributions are quite insignif- 
icant. It is clear that some states involve constructive 
and some destructive interference between the 1-impurity 
wavefunctions in different regions of space, but the spa- 
tial patterns are, not unexpectedly, considerably more 
intricate than the "hydrogen molecule" type states one 
might first imagine would form, with electrons living ci- 
ther directly between impurities or completely expelled 
from this region. This is of course due to the d-wave char- 
acter of the medium in which the quasiparticles propa- 
gate. For example, the LDOS is zero at the point halfway 
between the two impurities for the p-wave states, but it is 
quite small in the s-wave states as well. It is furthermore 
clear from the Figure that both s and p functions can 
have either constructive or destructive character, in the 
molecular sense. Note that the states are shown arranged 
vertically according to their eigenencrgies, but recall that 
the ordering of the s and p (T>i and T>2) states changes 
according to whether R is even or odd, as indicated in 
Figure |21 

In Figure we show the full energy variation of the 
LDOS spectra on sites near one of the impurities to illus- 
trate the expected widths of the resonances and the vari- 
ability with position. Note the surprisingly sharp high- 
energy resonances, which we discuss further below. On 
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the other hand, on certain sites in the configuration of 
FigureElsuch as (0,0), (2,2) and (4,4), virtually no LDOS 
weight at all is observed. In Figure H3 we show LDOS 
maps for several impurity separations along (110). Note 
that in the high-energy p states for R=(2,2) and (6,6), 
both of which are fairly narrow, nearly total destructive 
interference exists along the (110) direction outside of 
the two impurities, but there is substantial weight along 
the (110) direction relative to each of the two impuri- 
ties. As one moves the impurities apart, the tails of 
these wavefunctions in the (110) direction simply shift 
with the impurities. This suggests that these states may 
be somewhat narrowed because the quasiparticle is lo- 
calized between the two impurities, but only in the (110) 
direction; it may easily leak out along the (110) tails. 

2. R || (100) 

We noted above that for large separations the interefer- 
ence between quasiparticle wavefunctions vanishes much 
more rapidly with increasing separation in configurations 
with impurities aligned along a crystal axis R || (100). 
For small or intermediate spacings R < £o, however, the 
bound state splittings are just as large. Figure \7\ shows a 
spectrum for impurities separated by 6 lattice constants 
in the (100) direction. Although the low energy peaks 
are weak, there are nevertheless four well defined peaks 
as expected. The most striking feature of Fig. [3 is that 
the upper resonance is extremely sharp, far sharper in 
fact than a single impurity resonance at the same en- 
ergy! This is counter-intuitive based on our knowledge 
of the one-impurity problem: the T-matrix denominator 
S± defined in Eq. J5J) has an imaginary part proportional 
to the density of states of the clean d-wave superconduc- 
tor, so that the resonance width depends (approximately) 
linearly on the resonance energy. This is clearly not the 
case here. In the two-impurity problem where one impu- 
rity is at the origin and the second is at R, the T-matrix 
denominator is given by Eq. @, which can equivalcntly 
be written as 

V = det[vS - G (limp) (R, R,w)]detf (w), 

where T is the one-impurity T-matrix and G^ llmp ^ is the 
Green's function with one impurity at the origin. Thus 
sharp two- impurity resonances occur for exactly the same 
reason as in the one-impurity case, but because the one- 
impurity DOS at R is nonmonotonic in w, the resonance 
broadening is not necessarily proportional to the reso- 
nance energy. 

It seems intuitively clear that, because of destructive 
interference along the nodal directions, 2-impurity bound 
states could be formed in which quasiparticlcs are quite 
effectively trapped because they will be prevented from 
escaping via the long (110) tails of the individual impu- 
rity wave functions. In Figure [8] the spatial structure 
of the high-energy states confirms our intuition, since 
quasiparticlcs appear to be confined primarily to the axis 
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FIG. 7: Comparison of LDOS spectra p(r, uj) for various 
r with two impurities at positions (—3, 0) and (3, 0) (R = 
(6,0)). Tight binding band, A = 0.1, V = 10, fj, = 0. 



joining the impurities. The 45° tails of these 2-impurity 
states are supressed in all directions, and as expected the 
spectral features are even narrower than for the R || (110) 
configurations of Figure [SJ In the R = (2,0) case, the 
quasiparticle is nearly completely confined to the site be- 
tween the two impurities, while the wave function spreads 
out somewhat in the R = (6, 0) case. For these cases 
Gi(R, lo) vanishes as in the (110) case, so the wave func- 
tions have the same structure as in Eq. (|25I26|I . We 
note that when the peaks are sharp, the resonance state 
is essentially localized, that is the contribution to the 
wavefunction at resonance from the continuum has van- 
ishing weight. In general, one might have expected these 
states to decay as 1/r along the (110) direction and expo- 
nentially along the (100) direction, but they may decay 
more rapidly from interference effects. 

Thus far we have considered only a simple tight- 
binding band at half-filling, and one might worry that 
the existence of quasilocalized states was a consequence 
of the perfect nesting of this special electronic structure, 
and that such features are unlikely to be observed in 
real systems. We show below that this is not the case, 
and argue that a commensuration of dominant scattering 
wavevectors at the bound state energy is the important 
quantity, and that finding such a state depends on band 
structure and other details. 



VI. REALISTIC BANDS 

Until now we have focussed on the rather artificial 
symmetric tight-binding band case, both for calculational 
simplicity and to illuminate the unusual density of states 
phenomena driven by nesting features of the Fermi sur- 
face. Some, but clearly not all of these phenomena will 
survive away for a realistic band with particlc-hole asym- 
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FIG. 8: LDOS maps at resonant energies for R |] (100). Tight 
binding band, Ao = 0.1, Vo = 10, fi = 0. 
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FIG. 10: LDOS spectra for various sites as indicated. Upper 
panel for impurities at (-3,-3) and (3,3) (R = (6,6)), lower 
panel for impurities at (-1,0) and (1,0) (R = (2,0). Parame- 
ters for "realistic" band (see text), energies uj in units of ti, 
with Vo = 5.3. Note peak at lo = —0.18 in lower panel is gap 
edge, not impurity resonance. 
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FIG. 9: Solid line: BSCCO total density of states after Ref. 
[33l | vs w/ti. Dashed line: with d-wave superconducting gap 
of magnitude Ao = O.lti. Insert: Fermi surface at optimal 
doping. 



with t ...t 5 = 0.879,-1,0.275,-0.087,-0.1876,0.086 
and \ti \ = 0.1488eV". The density of states of this band 
in both normal and superconducting states is shown in 
Figure ED 



metry. Because of the intense current interest on STM 
studies of the cuprates, we now focus exclusively on a 
"realistic" representation of the electronic structure of 
BSCCO-2212, which we parametrize by adopting the 
tight-binding coefficients of Norman et al. |3^| : 

e(k) = t + 2ti [cos(k x ) + cos(fc. y )] + At 2 cos(k x ) cos(k y ) 
+2f 3 [cos(2fc x ) + cos(2fc y )] 
+2t4[cos(2k x ) cos(fcy) + cos^) cos(2fcy)] 
+At 5 cos(2fca;) cos(2fcy) (27) 



A. Direct spectra 

We first reconsider some of the impurity configurations 
we treated in Section to see how the LDOS patterns 
and bound states are affected by the electronic structure. 
The upper panel in Figure ED shows the distribtion of 
spectral weight on nearby sites for a configuration with 
impurities at separation R = (6,6). Comparison with 
the half-filled tight-binding model case shown in Figure 
[3] shows little qualitative change on nearby sites with re- 
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FIG. 11: LDOS maps for realistic band (see text), Vo = 5.3ti. 



spect to either peak positions, spectral weight or spatial 
distribution. On the other hand, comparison of LDOS 
maps in Figures |H1 and ^2 show that the longer-range 
character of the hole-like resonances has changed con- 
siderably. In other cases, particularly for small R, the 
changes are much more drastic. In particular, one does 
not generally see all four resonances, as also found by 
Morr and Stavropoulos[l8j . This is particularly true for 
resonances with R || (100), as illustrated in the lower 
panel of Fig. 1101 where the low energy resonances have 
completely disappeared on all sites investigated. 

To give an impression of the systematics of impurity 
resonance dependence on separation in the realistic sys- 
tem, and the configurations in which certain resonances 
are overdamped, in Figure^l we present a series of plots 
of LDOS spectra on nearest neighbor sites along the line 
joining the impurities for a range of separations with 
R||(100) and (110) directions. The number of peaks in 
each spectrum is variable, and for the closer separations 
normally only two peaks are observed (no peaks are ob- 
served for separation l!). As the impurities are moved 
farther apart, hybridization weakens and it becomes eas- 
ier to observe the full complement of four resonances. It 
is also clear that the bound state splittings for the real- 
istic band are decaying faster in the (100) direction, but 
that this splitting has not disappeared even for separa- 
tions as large as R = 13. This suggests that even in 
relatively dilute impurity systems the assumption of iso- 
lated impurities used to analyze recent STM experiments 
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FIG. 12: LDOS spectra for realistic band and Vo = 5.3ti on 
nearest neighbor site. Upper panel: impurities at (— R/2,0) 
and (.R/2,0) ( R = (R,0)), spectra taken at r = (R/2,1). 
Lower panel: impurities at (—R/2, —R/2) and (R/2, R/2) ( 
R = (R, R)), spectra taken at r = (R/2, R/2 + 1). 



may need to be reexamined. 

We note further that extremely sharp states occur fre- 
quently, for both even and odd separations; clearly the 
commcnsurability condition depends sensitively on the 
details of the band structure. Contrary to the results 
of Ref. 01 1 we nn( i that the state farthest from the 
Fermi level is frequently sharpest. In general, however, 
the spectra found here for the band of Ref. [33| are quite 
similar to those obtained in Ref. |l8j when direct compar- 
isons are possible. 

The condition for a true bound state [see Eq. ©] is 
satisfied at frequency u by G°(R, u>)f G°(R, uj)T = 1. 
Since this must be satisfied independently for real and 
imaginary parts of G (R, a/), sharp resonances only ap- 
pear for selected impurity separations. The product 
G (R,w)fG (R,w)f in Eq. @ is equivalent^ written 
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FIG. 13: Fermi surface of BSCCO - 2212 with constant 
energy surfaces at uj — 0.04 shown as small filled ellipses at 
the nodal points, qi, q2, q3 are wave vectors for which the 
joint density of states is large. 



e iq R G(k, w)T(w)(2(k + q, u)f(u). (28) 

k,q 

It was argued by Hoffman et al.|l3j|. in an analysis 
of disordered BSCCO samples, that the characteristic 
wavevectors found in the spatial Fourier transform of the 
LDOS are determined by peaks in the joint density of 
states X)kl m Gn(k, w)Im Gn(k + q, u>). It is interest- 
ing to ask whether these same (/-vectors are seen in Eq. 

At the small energies considered here, 3 distinct q- 
vectors contribute to the joint DOS [a fourth, q = 
does not produce any oscillatory R-dependence in Eq. 
(1281) ] . as illustrated in Fig. [ED (we neglect, for states suf- 
ficiently close to unitarity, the distinction between the 
tips of the quasiparticle constant energy contours and the 
nodal wave vectors). As a simple example, consider the 
realistic band model in the (100) direction, shown in the 
upper panel of Figure IT2l The R-dependence of Eq. I|28|) 
is straightforward since qi • (1,0) = q2 • (1,0) « 2.28/a, 
and q3 ■ R = 0, where a is the lattice constant. Naively, 
the standing-wave condition for a particle trapped be- 
tween the two impurities is (qi • R + 2?7o) = mt, where 
i]a = tan _1 (7rA r oVb) is the one- impurity scattering phase 
shift. On the other hand, for R in the (110) direction, 
we must simultaneously satisfy the commensurability re- 
quirements (qi • R + 2 770) = (q3 • R + 2i]o) = n7r i an d 
(q2 • R + 2t7q) = mir to form a resonant state, and we see 
that — as we observe in Fig. 1121 — sharp resonances occur 
much less frequently in the (110) than in the (100) di- 
rection. Quantitatively, the criterion for standing wave 
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FIG. 14: LDOS map for 1 strong repulsive impurity (Vb = 
5.3) at hole type resonance Qq = O.Ollti in system with 
"realistic band" (see text) a) without and b) with filter. 



formation is approximately satisfied for R = (3, 0), (7, 0), 
(11,0), (14,0) and R = (3,3), (11,11), which generally 
agrees with Fig. ^| 

We stress, however, that the relative success of this 
naive picture at making quantitative predictions is a bit 
surprising. There is nothing in our consideration to ac- 
count for the w-dependence of the T-matrix, or for the 
Nambu structure of the Green's functions. Furthermore, 
the approximation of the integral over q in Eq. I|28(l by 
a sum over a few dominant wavevectors is not expected 
to by justified at a quantitative level. Nonetheless, our 
considerations seem to indicate that the long-lived two- 
impurity bound states are derived from a few selected 
wavevectors. 



B. Filtering effects 

The discrepancies between the simple picture of a Zn 
impurity as a strong potential scatterer in a d-wave su- 
perconducting host and the LDOS measured near Zn 
impurities in BSCCO-2212 sam ples have been alluded 
to above. Martin and Balatsky[33| proposed that this 
problem could be resolved by noting that electrons must 
first tunnel through the BiO layer before reaching the 
Cu02 plane; applying the appropriate matrix elements 
for this process led them to a picture of a "filtered" DOS 
in which, in the simplest version, the STM tip samples 
not the LDOS corresponding to the atom directly under 
it, but rather to a sum of the LDOS on the surrounding 
4 nearest neighbor sites. With this ansatz, the LDOS 
pattern surrounding a Zn atom becomes, at a resonant 
energy of -1.5meV, rather similar to the experimentally 
observed one, with a bright spot at the center of the pat- 
tern, see Figure ITU 

We now point out the rather obvious fact that this fil- 
tering mechanism is characteristic of the presence of the 
BiO layer in the BSCCO-2212 system, and should there- 
fore be present in any STM measurement. If two nearby 
impurities are located via their resonant signals in such a 
measurement, the same filter should be applied to extract 
the true LDOS of the superconducting Cu02 layer. If it 
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FIG. 15: Comparison of true and filtered LDOS for two 



impurities at (-1,0) and (1,0) R = 
0.026ti, "realistic band" (see text). 



(2,0). V = 5.3, |fif| 



is found that the filtering mechanism works only in the 
case of isolated impurities, but the observed pattern in 
the 2-impurity case is quite different than that predicted 
by the simple filtered potential model, it must be aban- 
doned and more sophisticated explanations sought. For 
example, it will be interesting to pursue the alternative 
"Kondo" explanation of Polkovnikov et al.0] in the case 
of the 2-impurity problem. If one takes this model seri- 
ously, bringing two impurities close together should in- 
duce an RKKY interaction between the local moments on 
each impurity site, suprcssing the local Kondo screening 
and thereby weakening each impurity's scattering phase 
shift. One might then naively expect, in such a scenario, 
that bound state energies would generally be found at 
higher energies than in the isolated impurity case. Of 
course, one's intuition based on the 2-Kondo impurity 
problem in a normal metal is to be distrusted in this case, 
where the linear bare density of states already makes the 
1-impurity problem critical 
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VII. CONCLUSIONS 

In this paper we have explored a number of aspects 
of the quantum interference of impurity bound states in 
e?-wavc superconductors. We gave the exact form of the 
2-impurity t-matrix for two potentials separated by R, 
and showed that in general it has four resonances at fre- 
quencies ±f2i and ±fia which depend on R. In sim- 
ple situations, the eigenfunctions of the 2-impurity res- 
onant states can be constructed explicitly in terms of 
the eigenfunctions of the 1-impurity problem. Depend- 
ing on the impurity configuration and electronic struc- 
ture, some of the resonances are overdamped on specific 



sites or indeed sometimes over the entire lattice, leading 
to a smaller number of visible resonances in special situ- 
ations. On the other hand, in other situations resonant 
states were observed in the 2-impurity problem which 
were much sharper than their 1-impurity counterparts, 
in some cases occuring quite far from the Fermi level, 
contradicting one's intuition that these states should be 
more strongly damped. We have interepretcd these states 
as impurity "traps" in which quasiparticles are hindered 
by quantum interference, over surprisingly long lifetimes, 
from leaking out of the region between the two impuri- 
ties. 

The splitting of the bound state energies relative to 
the 1-impurity case was studied, and it was shown that 
the parity and energy of the 2-impurity eigenfunctions 
oscillate as a function of impurity separation. At asymp- 
totically large distances the splittings were shown to vary 
as ~ 1/R for impurities aligned along the (110) direc- 
tion and - exp(-R)/VR along (100). Systematic STM 
measurements of these splittings for isolated pairs of im- 
purities at different R were shown to provide a direct 
measurement of the spatial dependence of the Green's 
functions of the bulk superconductor. 

Finally we calculated the local density of states for 
2 impurities in a realistic band characteristic of the 
BSCCO-2212 system on which most STM experiments 
have been performed. The one qualitative difference rel- 
ative to the particle-hole symmetric case we examined 
earlier was the ovcrdamping of some bound states on 
any lattice site we studied for certain impurity config- 
urations. This makes it clear that even some qualitative 
features of LDOS spectra with two impurities will de- 
pend on details of the system in question. To extract 
information from STM when 2-impurity configurations 
are isolated will therefore require a careful fit to the- 
ory. We have made predictions for several concrete situ- 
ations which can be tested if such configurations can be 
found. In particular, we have calculated the density of 
states for realistic parameters corresponding to a Zn im- 
purity in BSCCO, and given results for both the direct 
LDOS and for the "filtered" LDOS proposed by Martin 
and Balatskv|32j to explain discrepancies in the standard 
model of Zn as a potential scatterer when compared with 
experiment. If the "filter" works for isolated single Zn im- 
purities but not for pairs of Zn atoms, it would be strong 
evidence in favor of an explanation for the Zn results in 
terms of residual induced local magnetism of the defect. 

We close by remarking that the solution of the 2- 
impurity problem may have important implications for 
the disordered N-impurity d-wavc superconductor and 
the interpretation of STM experiments. In particular, 
we have shown on the one hand that pairs of impuri- 
ties can give rise to trapped states which have great deal 
of spectral weight; on the other hand, interference from 
other impurities can destroy the characteristic pattern 
expected for an "isolated" impurity even when they are 
widely separated. A more thorough investigation of these 
questions requires a careful comparison with the many- 
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impurity system. One hint of the importance of the 2- 
impurity states in the fully disordered system comes from 
the study of the perfectly nested band, where we find that 
many of the unusu al sy mmetry-based features of the to- 
tal density of states [33 are reflected already in the simple 
2-impurity problem as well. In addition, we have investi- 
gated the effects of self-consistent treatment of the order 
parameter on the results above, which were all produced 
assuming homogeneous Ak- We find that, although spec- 
tral weight is shifted by the order parameter supression 



around the impurity site, in general away from the Fermi 
level |39j. the LDOS patterns are rather weakly affected. 
We will report in detail on these findings elsewhere |3^l|. 
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